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Abstract. We consider a mixed stochastic differential equation involving 
both standard Brownian motion and fractional Brownian motion with Hurst 
parameter H > 1/2. The mean-square rate of convergence of Eulcr approxi- 
mations of solution to this equation is obtained. 



Introduction 

The main object of this paper is the following mixed stochastic differential equa- 
tion involving independent Wiener process B and fractional Brownian motion 
with Hurst index H £ (1/2, 1): 

(1) Xt = x„+ f a{.s,x,)ds+ f b{s,x,)dWs+ f c{x,)dBf, te[o,r], 

Jo Jo Jo 

where the integral w.r.t. Wiener process is the standard Ito integral, and the inte- 
gral w.r.t. fBm is the forward stochastic integral. The questions of existence and 
uniqueness of solution for equations of such type were considered in [7l [9l [5l [11] . 

Such mixed equations arise in different applied areas. In financial mathematics, 
for example, it is often natural to assume that the underlying random noise consists 
of two parts: a "fundamental" part, describing the economical background for a 
stock price, and a "trading" part, coming from the randomness inherent for the 
stock market. In this case the fundamental part of the noise should have a long 
memory, while the second part is likely to be a white noise. 

Due to a wide area of applications of equation ([T]), it is important to consider 
certain numerical methods to solve it. We use here the most popular and probably 
the simplest method of Euler approximations: one takes a uniform partition of 
the interval, where the equation is being solved, and replaces differentials by a 
correspondent finite differences. There is a vast literature dedicated to numerical 
methods for stochastic differential equations driven by the Wiener process, we refer 
to classical monographs 8 and j6] for an overview of the subject. There are also 
several papers dealing with discrete time approximations for stochastic differential 
equations with fractional Brownian motion, for example, [TUJ [T^l [3] . 

The main difficulty when considering equation ([1]) lies in the fact that the ma- 
chinery behind the two stochastic integrals is very different. The Ito integral is 
treated usually in a mean square sense, while the integral with respect to fractional 
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Brownian motion is understood and controlled in a pathwise sense. The mixture 
of two integrals makes things a lot harder, forcing us to consider very smooth coef- 
ficients and to make delicate estimates. 

The paper is organized as follows. In Section 1, we give basic fact about for- 
ward and Skorokhod integration with respect to fractional Brownian motion and 
formulate main hypotheses. In Section 2, we define Euler approximations of ([T]) 
and establish some uniform integrability results for them. Section 3 contains the 
main result about rate of convergence of Euler approximations for equation ([l}. 
Unsurprisingly, the rate of convergence appears to be equal to the worst of the 
rates for corresponding "pure" equations, i.e. the mean-square distance between 
true and approximate solutions is of order 5^1"^ V 5^^^^ , where b is the mesh of the 
partition. 

1. Preliminaries 

1.1. Fractional Brownian motion and stochastic integration. In this section 
we give basic facts about the stochastic calculus for fractional Brownian motion. A 
more extensive exposition can be found e.g. in [4l[l]. 

Fractional Brownian motion (fBm) is by definition a centered Gaussian pro- 
cess with the covariance 

E[i?fi?f]=i(t^^ + ,s2«-|t-.r),t,.>0. 

It has a version with almost surely K-Holder continuous paths for any k < H. 
For H e (1/2,1) (the case we consider here) it exhibits a property of long-range 
dependence. 

Let -/j|^[0, T] be the completion of the space of continuous functions with respect 
to the scalar product 

{L9)h= II f{t)g{s)i^{t,s)dsdt, 

J J[O.T]^ 

where ip{t, s) = H{2H - 1) \t-s\^"^^. Denote also by ||/||^ = ^(/,/)^ the 
corresponding norm. 

Now we recall the notion of stochastic derivative. Let infinitely differentiable 
function F : M" ^ R be bounded along with derivatives. For a smooth functional 
G — F{B^, . . . , B^), where ti, . . . ,tn G [0, T] the stochastic derivative is defined 
as 

n 

i?,G = ^F;jB,^,...,i?/^)I[o,t,](s). 

fc=i 

The Sobolev space B^'^ is the closure of the space of smooth functionals with respect 
to the norm 

\\G\\l,^E[G']+E 

The Skorokhod, or divergence, stochastic integral is the adjoint to the stochastic 
derivative in the following sense. Let the domain dom(5 of the divergence integral 
be the space of random processes u g L2{^, ^2^(0, T]) such that 

E[{DG,u)^]<Cu\\G\\^^^^^ 



\\DG\\l 



RATE OF CONVERGENCE OF EULER APPROXIMATIONS OF MIXED SDE 



3 



for all G e E)^'^. Then the divergence integral 

Jo 

is defined as the unique element of L2{i~l) such that 
(2) E[{DG,u)^] = E[GSiu)] 

for all G e D^'^. It is worth to remark that dom(5 contains the space D^'^(L^[0, T]) 
of processes such that 



I"IIh;1,2 = E 



ll^llff + / /// E[DtUyDsUz]ip{t, s)ip{v, z)ds dt dz dv 



is finite. Moreover, for such processes 

(3) E [Siuf] < \\u\\l.^,^ 



2 



The forward integral with respect to fBm is defined as the uniform limit in 
probability 



u,dB^ = lim / ^+'^ 
6^0 .L £ 



dB" = lim / u,^^ -ds. 



provided this limit exists. It is well-known (see e.g. [T]) that if for u £ D^'^(L|^[0, T]) 

\DsUt \ s)dsdt < oo, 



'[0,T]2 

then the forward integral exists and is equal to 

(4) / utdB^ ^ [ utSB^ + 11 DsUt'ip{t,s)dsdt. 

Jo Jo JJ[o,TY 

1.2. Assumptions. The following hypotheses on the ingredients of equation ([T]) 
will be assumed throughout the paper. 

(A) The functions a and 6 are bounded together with their derivatives o^, h'^: 

Ht,x)\ + \h[t,x)\ + \a'^{t,x)\ + \b',[t,x)\ < K- 

(B) the functions a and b arc uniformly {2H — l)-Holder continuous in time: 

\a{t, x) — a{s, x)\ + \b[t, x) — b[s, x)\ < K \t — s|^"^ ; 

(C) the coefficient c is bounded together with its first and second derivatives 
and uniformly positive: 

< c{x) + c{x)-'^ + \c'{x)\ + \c"{x)\ < K. 

Here if is a constant independent oi x, t and s; 

(D) the Wiener process W and the fractional Brownian motion B^ are inde- 
pendent. 

In what follows C will denote a generic constant, whose value might change from 
line to line. To emphasize dependence on some variables, we will put them into 
subscript. For a random process X we denote its increments by Xt^s ~ Xt — Xs- 
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2. EULER APPROXIMATIONS AND AUXILIARY RESULTS 

For iV > 1 consider the following partition of the fixed interval [0, T] : {0 — i^o < 
vi<--- <VN = T, 6 = T/N}, Vk = kS. 

The Eulcr approximation for equation ([T]) is defined recursively as 

where ISWk = W^^^^,,y^, = ^^k+i '^k' "^^^ initial value of approximations is 

Set = max{n : < u}, = I'n^ , and define continuous interpolation by 
Xt = Xf. + a{ti,X^sJiu - ti) + + c(Xf.)s5, , 

or, in the integral form, 

pU nU nU 

(5) Xt=Xo+ a{tixfs)ds+ b{ti, Xf,)dWs + c{Xfs)dB^ . 

Jo "Jo " Jo " 

The following lemma is a discrete analogue of the Gronwall inequality. 
Lemma 2.1. If a non-negative sequence {xn,n > 1} satisfies 

Xn+l < Xn{l +K5)+K5. 

Then 

Xn < (xo + l)e^^". 
The following two lemmas are technical. 
Lemma 2.2. For s < Vn, n>l, one has 

n-l 

(6) = c{X^s) n (1 + <i''k,Xtj5 + b',{i^k,X'jAWk + c'iX'jAB^) 

(the product is set to 1 when the upper limit is smaller than the lower). 

Proof. Clearly, DgXl^ = if Vn < s. Now observe that DgAB^ = I„=„5. Hence, 
for n — n^g, we have 

Further, for n > nf we can write 

D,X^^^ = Ds{xt^_^+a{yn^uXt^^_;)5 
+ 5K_i,Xf_JAT4^„_i +c(Xf_JAS,f_i) = (l + a^K^i, J<5 
+ 6;K_i,x5_JATy„_i + c'{Xt^_^)AB^_^)D,Xt^_^, 
and deduce © by induction. □ 
Lemma 2.3. For any M > it holds 

{N~l 
M J2 ii^Wk f + (ABf )2) \ < Cm 
k=0 

for all N large enough with Cm independent of N . 
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Proof. Using independence of W and , we then can write 



exp 



( N-l 



fc=0 



N-l 



exp <^ M i^Wk 



k=0 



'N-l 



Mexp ^(Ai?f)2 



fc=0 



N~l 



N-l 



l/N 



< n E [exp{M(AW^fe)2}] ( J] E [exp {MN{AB^)^}] 

k=0 k=Q 

= C(l - 2Af(5)-^/2(^ _ IMNS^^Y^I'^ 
= C(l - 2MT/iV)-^/2(^ _ 2AfT2^iVi-2//)-i/2^ 

where the last equahties hold provided 2MT/N < 1 and 2MT^^ N^^^" < 1, which 
is true for all N large enough. Observing that 

we get the desired boundedness. □ 

Now we are ready to prove that the moments of Euler approximations as well as 
of their stochastic derivatives are uniformly bounded. 

Lemma 2.4. For any p > one has 

(7) E[\D,XtJ]<Cp 
for all s e [0,T], n < N, with Cp independent of 6. 

Proof. It is easy to see from ([6|) that the left-hand side of ([7]) is finite. Therefore, 
it suffices to establish boundedness only for N large enough. 
Introduce the following notation: 

Qk = \l + ak5 + bkAWk + CkAB^\,^u - |AT4^,| + |ABf | , 
dk^a{yk,Xt^)5 + b[yk.XtjAWk.Ak = dk + c{Xt^)AB^ ^Xt^^^^^^. 

Fix a small positive constant 7 (its value will be specified later to satisfy our 
needs). Put A = {ik < 7}- 

E [|i?sXir] = E [\DsX'J''Ia] + E [\D,XtJln\A] ■ 

Step 1. First we estimate E [jD^X;^^ ■ Write 



E[\D,XtJtn\A]=Y.^ 



ixf^r n Qfc^7.<7 n Qfci 

7fc>7 

\k</B keB / 



where the outer sum is taken over all non-empty B C {nf, -I- 1, . . . , 
Observe that 

\akS + bkAWk + CkABl^\ <K{S + jk), 
so this expression does not exceed 1 whenever S < 1/{2K) and 7fe < 7 < l/(2i^ 
and we can write 

6fel7,<7 < exp {ak6 + bkAWk + CkABj^} . 
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For 7fe > 7 we estimate simply 9^ < exp {K{S + 7fe)}, therefore, 



< 



E[\D,X',Jln\A] 
expJ ^K(5 + 6feAM^fc + CfeASf) I JJ Qk^^.y^ 

\ [k^B ) keB / 



<CpE exp\pJ2iKS + bkAWk + K\ABj^\)\Y[ 

[ kiB ) keB 



fk>l 



<Cp E 



exp I 3p J2 bkAWk 

{ k^B 



expispi^ 5] |ASf|| 



xE 



ik>i 



.keB 



1/3 



By the standard properties of the stochastic integral with respect to W, 



(8) 



exphpJ2bkAWk \ 

I kiB ) 



E exp <^ ^ V52<5/2 
< exp {5p^K'^NS} = exp{5p'^K^T} . 



Now estimate, using the Holder inequality. 



< E 



N-l 



(9) 



l/N 



exphpJ^CkAB^ 

{ k^B 

CN-1 
fc=0 



exp\3pKj2 |A^f I 



fc=0 



< E 



3pKN\AB^\ 



Further, 



■Q e^pK-ykj 



7fc>7 



LkeB 



< E 



< E 



\6pK^ 2 
exp < 2^ 7fc 

I k=a 



jeB 



1/2 



jeB 



1/2 



where the last inequality hold for all N large enough thanks to Lemma 12.31 To 
estimate the last expectation, recall that W and are independent and take first 
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the expectation with respect to W: 



n ^7j>7 

jeB 



< E 



[]2<i>((|Ai3f |-7)<5-i/2 

jeB 



(10) 



<2"(s) I [|e $((|ABf|-7)r 



1/2 



<i>((|A<|-7)<5 



-1/2 



«(B)' 



l/n{B) 



where rii^B) is the number of elements of B. We spht the inner expectation 
into parts where \AB^\ < j/2 and \AB^\ > 7/2. For \AB^\ < 7/2 it holds 
$((|AB^| -7)(5-i/2) < $(-7(5-i/2/2)"(^), also we have P{\AB^\ > 7/2) < 
2<i>{-jS-"/2), hence 

E h{{\AB^\ - 7)<5-i/2)"^''^] < $(-7,5-i/2/2)"(B) + 2^~-fd-" /2) 

< g-7'<5"'n(S)/8 ^ g-7'<5-'"/2 <- g-C-,«(B)Ar _^ g-C^JV^" ^ 

Plugging this into (fTO)) . we get 





< E 







n2<i>((|Ai3f|-7)'5-^/') 

jGB 

and combining this with ([8]) and (|9]), we arrive to 

E[\D,X',Jln\A] 



\ S B J 



For TV large enough it holds C(C^iV — 1) > C^iV (naturally, with different constants 
in the left-hand and in the right-hand sides), so the first sum is bounded from 
above by 



< (1 + e-^-'^'f - 1 = exp {TVlog (1 + e-*^-^) } - 1 



< exp [CNe-^''^] - 1 < exp [C^e-^'-^] -l<C^t 
Similarly, the second sum is bounded by 



-C^N 



^CN-C-,N^ 



Since H G (1/2. 1), this implies 

E[\DsXtJln\A\ <Cp, 



7^ 
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for all TV large enough. This expression vanishes as — >■ oo, hence we get the 
boundedness of E [|Z3sX^^ ■ 

Step 2. Now we turn to E [\DsXI^ \^1a] ■ If we take 7 < K''^ /'i and 5 < K'^/S, 
then|Afe| < 2K-'^/3 on A and \c{xl^^ J- c{X^J\ <2K~^/3. But c{X^J > K-\ 
so c{X^^^_^)/c{Xf^J G (1/3, 5/3), which allows us to write by the Taylor formula 



(11) 



log 



Ck 



where \R'f.\ < CAl. Similarly, on A 

log Qk - log (1 + akS + bkAWk + CfeASf ) = akS + bkAWk + CfcASf + R'^ 

with \R'i^\ < CA^. Plugging into this formula the expression for CkABj^ from (fTTj) . 
we get 

c(^^ ) 

= ,"s^\ exp {akS + PkAWk + Rk} , 

where Rk = R'l - R^ «fc - - Cka{XtJ/c{XtJ, /3fc = hk - Ckh{XtJ / c{Xt^) . 
Now we can estimate 



£[\DsXIJ1a] =E 



n-l 



,{xtsY n 0'!^ 



k=ni 



z{Xt^_Jex^lpY, {akS + PkAWk + Rk) 



k=ni 



< a, 



n-l 



exp <^ 2p PkAWk 

\^ k=nl 
( n-l 

exp <^ 2p' J2 



N-l 



k=nl 



exp < 2p ^ \Rk 

I /c=0 



1/2 



exp 



1/2 



< CpE 



exp 



N-l \- 
k=0 ) . 



where the last holds for all N large enough due to Lemma [2731 This completes the 
proof. □ 

Lemma 2.5. For any p > one has 

(12) E[|X/n<Cp 

for all t € [0, T] . 

Moreover, 

(13) 

for any t G [0, T], 



- A 5 



< CpSP/'^ 
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Proof. It is enough to prove this for p = 2m, m G N. We first prove (I12|) for t = Vn, 
using an induction by m. 
Start with m = I. 

Denote a„ = a(i/„, X^^), 6„ = 6(t'„, X^^), c„ — c(X^^) and write for 6 G (0, 1/2) 
by Jensen's inequahty 



< (E 



5 n2 



E [(6„AT4^„)2] + E [(c„AS,f )2l 



+2E [X5^6„A1^„] + 2E [6„c„ AW^„ ASf ] + 2E [xf^ASf ] )e2^ + C5 



< E 



< E 



CS + C6'^" + 2E [xf,^AB^] )e^^ + C6 



X. 



5 n2 



By ^ and (O, we can write 



E [< AB,f ] = an /" E [t - sf"-^dtds 

Jo Ju„ 

<C f "^\t- iynf"-^dt < CS. 



Then by Lemma 



X 



S 



< X^e^*" < Ce^'^ < C, 



as required. 

Now let m > 2 and for I < m 



< C21. 



In the further estimates constants may depend on to, but not on n. 
Observe that by the Jensen inequahty for S < I 



{a + bf"' < (1 - (5)l"2mQ2m _^ ^l-2m^2m 



whence 
(14) 

therefore 



[a + 6)^" < a'"(l + C^5) + C^b^"'S 



2m rl — 2m 



< E 



2m 



?l-2m 



< E 



«^ + 6„AT4^„ + c„Ai?f (1 + CS) 
{X^^^ + b„AWn + c„ABf (1 + CS) + CS. 
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Expand the power in the first term and consider a generic term of this expansion 
(without a coefficient): 



X, 



X6„AT4^„)'(c„Ai3,f)'= 



= E 
= E 

Thus, we can write 



2V2(z/2) 



X „ + hn/XWn + C„AB, 



H 



2m 



2m m 



2m 

fc, 2j, 2m- k- 2j 



X 



g \ 2m — 2j — k 



6?/c^(Ai?f)* 



20 j\ 



5' 



where 



a + 6 + c 
a, fe, c 



[a + b + c)\ 



a\b\c\ 

is a trinomial coefficient. 

For k — 0, j > 1, the terms of this sum are bounded by CS by the induction 
hypothesis and boundedncss of 6„,c„. 

Further, for fc > 2 



< CE 



X 



S ,2m-2j-k 



AB, 



HI 



< CE 



g \ 2m 



l/A 



AS, 



<c(l + E[(xfj' 



6^" <C{1 + E 



1/t, 

(5 \ 2m 



where A — 2m/ {2m — 2j — k),rj = A/(A — 1); here we have used an obvious estimate 



(15) 



2m 



l/A 



<1 + E[«)]' 



Now estimate the term with k — \, j — Q, using formula ^ 
'c„AB,f = / / E D.AiX 



^i/. 



\2m-l 







r 


L ' 


Jo 





c{X^J] %l^{t,s)dtds 

5 



< C 1 



((2m - l){xtS"'-'c{Xtj + «)^"-^c'«)) 
xH{2H ~ l)(t - s f"-^dtds 

{t~:y,,)^"~'dt<c{l + E[{x'^S 



g \ 2m 



s. 



Here, as above we have used the Holder inequality, inequality p5|) and boundedness 
of moments of the stochastic derivative. The terms with A: = 1 , j > 1 are estimated 
similarly. 

Collecting the estimates, we get 



X 



< E 



{Xlj^"" {l + CS) + CS, 
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ic B" r 



SO by Lemma [2. II we get the desired boundedness. 
Now write for s G [vn, i^n+i) 

E - l"] < Cp(E [\anis - ur,)n + E [|6„M^,,,jn + E 

which gives ([13]) and together with ([12]) for t = imphes (IT2|) for alH G [0, T]. □ 

3. Rate of convergence 

Now we are ready to prove the main result about the mean-square rate of con- 
vergence of Euler approximations. 

Theorem 3.1. Euler approximations ([5|) for the solution of equation Q satisfy 

E UXt-XfY] < C{S + S^"-^). 

Proof. Define ip{x) = c{z)^^dz. It is clear that 

K^^ \x ~ y\ < \i^{x) — ip{y)\ < K \x — y \ . 
Write by the Ito formula 



where 



Similarly, 



a{s,x) b{s,xYc'{x) b{s,x) 



c{x) 2c{xy 



c{x) 



where 



\xl)(^{a{s,Xl) - a{tlxf,))ds + {b[s,X'^) -b{tlxf,))dWs) 



+c~\xt){c{Xt) - c{Xts))dBf + - b{tlxt.)fds 

f (af ds + 6f dWs + cidBf + dl ds) . 
Jo 



Estimate 



< G 



atds 



<G f E [{alY] ds 
Jo 

c{x!)-^a{s,Xf,)^aitlxf,)y + {a{s, x!) - a{s, Xfs)f 



ds 



X 



S \2' 



ds] <CS'"-\ 



Similarly, 













Jo 



<G f E[{dif]ds<GS^"-\ 
Jo 
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and using the Ito isometry, 



J^bidWs^ = J*E[{bl)']ds<CS' 



2H-1 



Further, by dJ) 



cf SBi 



?f + // D,ciij{s,u)dsdu^:I'{6,t)+I"{S,t). 
By the chain rule for the stochastic derivative, 

-(c(Xf ) - c{X^s))D^x! + c-\xl){c!{Xl)D^Xl - c'{Xf,)D^Xf,) 



c{XtY 



(c(Xf ) - c{Xfs))DuX! + c-'{x!){c'ix!) - c'iXfs))D^X^ 



c{Xlf 

+c-\x!)c'{Xfs)DuXl,s) Di{u,s)+D2{u,s) + D3{u,s). 



Now 



E [Diiu,sf] < CE [(c(Xf) - c(X/.))'p„xf)2 
Similarly, 

E[D2{u,sf] <CS. 

Further, for u < 



< CE 



{x!^x^sy 



< cs. 



x' 



and 



DuXfs 



(aUtlxf,){s - tf) + b',{tixfs)W,^ti + c'iXf^jB^,. 



1/2 



xE 



E[D3{u,sf] <E\{DuX^^ 



1/2 



< CS; 



for u e [if, s) 

D3iu,s)^ciXfs) 
and D'i{u, s) = for u> s. Thus 

E [I"{5,tf] < C 



Di{u,sy + D2iu,sy 



+ Dz{u^ s)^I[o,t«] {u) tp{s, u)du ds + CE 



Jt'i 



-0(5, u)du ds 







{I>- 


s 1 







2H-1 



<C\6+[ \s~a ds] ] <Ci6 + 6^"-^). 
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By 111) 

E [I' {6, tf] < [ E [cf ci] V'(s, u)ds du 

J[0,t]2 

+ / E[DucI DgC^] 'ip{s,u)dudv ds dz. 

JlO,t]i 

The first term is estimated hy CS using that 
E[cfca<E[(4f]+E[(cff]<CE 



< CS. 



In the second, we write -Ducf = Di{u,v) + D2{u,v) + D^{u,v) and similarly for 
Dsci. For Di,D2 we use the Cauchy inequality and the above estimates to get 
a bound of C6. This also works for when u ^ [ti,v) and s ^ [ti,z). Two 
remaining terms for D3 are similar, take e.g. 

pt pt pt pV 

/ / / / E [Dtjcf i^scf] |s — u| k ^ ''^l dudvdsdz 
Jo Jo Jo Jtl 

<C5^"-^ [ f f E[\Dsci\]\z-vf"~^dvdsdz. 
Jo Jo Jo 

Again, if s ^ [tliz), the integral can be estimated by C5^^^; for s £ [tl,z) we get 
pH-i Ultimately, 

and adding this to the previous estimates, we get 

E[(G*)2] <C{d + S^"-^). 

So we can write 



(a(s, X,) - a(s, Xi)f + (/3(s, X,) - X^Y ds + E [[Clf] 



< C / E 



< C / E 

Jo 



By the Gronwall lemma. 



hence 



{Xt - Xf 



,1/2 



□ 



< 



E {Xt^xfy <cis + s^"-^), 

as required. 

Remark 3.1. The obtained estimate can also be written as (E 

C{S^^^ V S^^^-^), so the mean-square rate of convergence for the mixed equation is 
the worst of the two rates for "pure" stochastic differential equation with Brownian 
motion, C5^/^, and with fractional Brownian motion, C5^^^^. As long as these 
estimates for pure equations are sharp (see O [E]), we get that in our case the 
estimate is sharp as well. 
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An interesting observation is that the value of the Hurst index where the rate 
of convergence changes is i7 = 3/4. From 2 it is known that the measure induced 
by the mixture of Brownian motion and independent fractional Brownian motion is 
equivalent to the Wiener measure iff iJ > 3/4. So in this case it is perhaps natural 
to expect that the rate of convergence of Euler approximations is the same as for 
Brownian motion, and this is exactly what we see here. 
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